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ABSTRACT 

A new theory of thermal history after inflation is outlined by taking into account 
the effect of parametric resonance caused by inflaton oscillation. A typical outcome 
is existence of a high temperature phase, that makes possible GUT baryogenesis 
consistently with low enough gravitino abundance. 
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It is natural, and even inevitable as the fundamental microscopic theory advances, 
to ask how things in the universe are as they are. The classic example in cosmology 
is how the element abundance is as it is; the problem of nucleosynthesis. In recent 
times we went quite further along these lines: we now even talk about origin of 
matter-antimatter asymmetry, creation of the spacetime itself, and so on. What I 
would like to discuss in this talk is how the bulk of entropy in our universe, the 
universe at least in the observable portion to us, has been created. The problem 
might sound purely accademic, but actually has very practical aspects with regard 
to the problem of baryogenesis at the Planck or the GUT epoch, as explained later. 

We shall discuss the entropy generation in the paradigm of inflationary universe 
models. This puts our task to the logical extreme, since at the end of the exponential 
expansion of inflationary stage our portion of universe becomes radiation free: what is 
left from inflation is simply oscillation of the inflaton field (coherent and homogeneous 
field that causes inflation) around the minimum of inflaton potential. Thus from the 
entropy point there is essentially nothing at the start. We must understand how the 
entire entropy of our universe is generated originating from inflaton oscillation. This 
problem has been often called reheating problem, although it is somewhat misleading 
to put the word "re" because it does not matter whether there is an epoch of thermal 
equilibrium prior to inflation. The elementary picture of reheating until the middle 
of 1994 was as follows : the coherent and homogeneous inflaton oscillation can be 
regarded as an aggregate of particles at rest with a precisely tuned phase relation 
and these nonrelativistic particles independently decay according to the probability 
law given by the simple Born decay formula. Assuming these produced particles 
immediately interact strongly, one may hope to achieve rapid thermalization among 
created particles. 

Recent works on this subject made it clear that this elementary picture is grossly 
wrong if the initial oscillation amplitude of the inflaton field is large [0], [^, 
12 1 , [13 1, [^], [0. Particle production under a large amplitude and almost periodic 



oscillation is dictated by many instability bands only one of which gives the Born 
decay rate, a valid approximation under the small, but not under the large amplitude 
limit. Existence of infinitely many bands under a periodic coefficient function is 
known as the phenomenon of parametric resonance |]ll|, 0. Incorporation of this 
new source of many instability bands gives rise to much faster inflaton decay and 
much more copious particle production. Thus the thermal history right after inflation 
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is quite different from what was anticipated previously: there exists a prolonged 
epoch of much higher temperatures than expected from the Born formula 0. It 
opens up a new possibility of baryogenesis consistently with the gravitino problem, 
as discussed later. A different picture of the reheating has been advocated by Linde 



in the preceeding talk [TQ • 

As a model of inflation we take the simplest; chaotic inflation with a mass term 
potential | m| . At the Hubble rate equal to the inflaton mass, damped oscillation 
starts. To be consistent with COBE anisotropy of the microwave fluctuation we 
may take the inflaton mass in the range of = 10^^ GeV. The inflaton field has 
coupling to ordinary matter and gauge fields, which we generically denote by ip. 
Only bose fields are relevant to the following discussion. These couplings may be 
either of the quartic ^^v^^ or of the Yukawa C,'^'^ type. Ratio of these two terms is an 
important quantity to the rest of discussion, but without much guidance from the 
fundamental theory we take a simple supersymmetric model with a superpotential 
of the form, W = |m^^^ + "^fi'^V^^? which gives an ordinary potential for the C, 
coupling, = I g'^i'^v'^ + \ S'^^g^V'^ • Fixing the ratio of the quartic to the Yukawa 
coupling terms to roughly of order g^/m^ is important, but a precise ratio number 
given here does not matter to the rest of conclusions. 

For the time being, let us ignore time variation of the inflaton amplitude due to 



the cosmological expansion. We find it best to exploit the Schrodinger picture |13 



in order to clarify the nature of quantum state under the periodic perturbation of 
the inflaton field. In this formalism any quantum state can be expressed in terms of 
a direct product of the state vectors \iljit))k which obeys the Schrodinger equation 
of a periodic frequency ujl{t). It turns out [|13| that the quantum wave function is 
solved with a Gaussian ansatz; 

where the time dependent function Uk{t) obeys the same classical oscillator equation 
as the mode equation corresponding to the wave vector k; 



+ cu^ (t) Mfc = , uoi{t) = k' + g'e + gm^^ 



with ^ = ^0 cos[mg(t — to)]. Moreover, an initial choice of the quantum state can 
be translated to a choice of the boundary condition to this classical equation. We 
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took the ground state since inflation pushes out essentially all matter out of horizon. 
The initial vacuum at t = to then corresponds to u{t) ~ (c(j(to)/7r)^^/^ e*'^*^*"^ . 
An important point here is that there is no freedom to impose the boundedness at 
t = oo, quite unlike the Bloch wave solution in condensed matter physics (time is 
replaced by coordinate there). 

A mathematical theorem p[ tells that with periodicity of uj'^it) alone, and ir- 
respective of any details of this function, there are infinitely many bands of in- 
stability separated by stability bands in the two parameter space that character- 
izes co'^(t), k and the oscillation amplitude ^o- Instability grows exponentially with 
time: the asymptotic form of solution is u{t) e^*x (periodic function) with 
/i > 0. Furthermore, within each instability band denoted by n there is an eigen- 
mode along a curve in the parameter space, k = fc„(^o) ; ('^ = 1 > 2 , 3 , ■ ■ ■), that 
has periodic zeros at tj (i = 1 , 2 , 3 , ■ ■ ■), M„(ti) = 0. This eigenmode approaches 
k —>■ ^ as .^0 ^ 0. This remarkable fact has a profound implication on 
the behavior of the quantum state: the phase factor ^ ln|M(t)p in the wave 
function (g| ip{t) ) diverges at these tj. All off-diagonal density matrix elements 
Pil il' ii) = (<?! i'i't) ) (^'l i^it) )* 7 (97^ 1') 5 vanish, when time averaged over short 
time scales of order We take this gnal to indicate an almost classi- 

cal behavior in the quantum system under the periodic perturbation. The coarse 
grained density matrix thus defined is diagonal and positive definite: in the harmonic 
oscillator Fock base of \n) = u {n + ^)\n) , (n = , 1 , 2 , ■ ■ ■) , Pn,m = Pn ^nm , 
with 

00 00 

EPn=l, E(Pn)^<l- 

n=0 n=0 

Thus the system is effectively described by a classical probability distribution func- 
tion p^. 

This classical instability implies that the wave function of the quantum system 
has an exponentially growing Gaussian width (g^) oc |'u(t)P , which means that 
higher energy levels are excited with time, to be interpreted as particle production 
in the particle language. In this way one can precisely compute produced number 
of particles in terms of the density matrix element, resulting in (A^) = Pqq — 1 for 
each mode. The formalism manifestly respects unitarity, and one may devise useful 
approximation for many quantities. For instance, the particle number distribution 
P2n,2n is expressed in terms of the average number of particles (A^) alone, which as 
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n ^ oo, has an asymptotic form of 

P2n ,2n — p2n ,2n ~^ / ; /—— ■ 

in + 1) 

I would like to stress that this distribution is rather broad, as evidenced for instance 



by the calculated dispersion ^ {N'^) — {N)"^ / (N) ~ 1 . The pseudo-entropy defined 
by — tr In is calculated as ~ In(A^) + 0.09 using the asymptotic formula, which 
grows logarithmically with the degrees of freedom (A^) as usual. 

With cosmological expansion one has to consider damping of the inflaton am- 
plitude ^(t). For the Hubble time 1/H much larger than the oscillation time scale 
one may deal with averaged quantities over one oscillatin time. In the rest 
of discussion we take ^ to mean the time-averaged amplitude. In the model of our 
concern we may separate the initial and the final phase of inflaton oscillation, de- 
pending on the value of r/ = g^/m^ . At the very start of the ^ oscillation this value 
is of order rjo = g^o/ni£, ~ 10^(7, which is large unless g is extremely small. Both 
the large rj and the small rj regions are approximately described by the standard 
Mathieu equation; 

— + {h - 29 cos{2t))u = . 

For instance, in the large rj limit initially relevant, the quartic inflaton coupling 
dominates and h = {k'^ + = (7^^^/(4m|) = r = m^t . Thus 

the relevant instability band is limited to h — 26 {= k'^/m'^) > 0. Although much 
is known on instability of the Mathieu equation [Q, we had to improve much more 
in the large and the small 6 region, especially on analytic form of the growth rate 
X{h,6) in the late time behavior of e'^™'«*. Detailed results of this investigation are 
left to Ref 0, [0], but let me mention a few key points. 

First, in the large amplitude limit one may use analogy to scattering problem in 
quantum mechanics, to analyze the growth rate in instability bands. Starting at the 
initial "position", to, one considers infinitely many barriers of the same height, each 
of which can be replaced by an inverted harmonic oscillator in the large 6 limit. Only 
the growing component is amplified each time of the barrier crossing, and one may 
obtain an analytic formula for the growth rate by diagonalizing the transfer matrix 
over one barrier crossing. The result is summarized ^ as follows: with e = ^ — 1 , 

A = - ln(v^+ V:^^) , z = (l + e""^') cos^t/^, 

TT 



5 



^ = + v^e ln(7r 6^'^) + - $5 ln[ -—- 

2 ^ ' 2 ^r((l + zv^e)/2 

The instability island that appears as e varies with a fixed 6' is characterized by 
X > 1, separated by stability regions of a; < 1. The maximal rate A is ~ 0.28 at 
e = 0, and as e increases from 0, the rate decreases. The averaged value taken over 
a large region of 6' is ~ 0.15 for the first island of smallest e. Since it is difficult to 
deal with each instability band separately and then add them together, we use this 
average rate and multiply the phase space volume of 



0\ 



27r2 



contributing to this rate. The average energy of created ip particle is of order y gm^^ , 
which is initially very large ~ ^/g ■ 10^^ GeV . 

On the other hand, in the the small amplitude limit the Yukawa coupling term 
dominates and the Mathieu parameter becomes h = Ak'^/m'^ , 9 = 2gC,/m^ = 2r] . The 
decay of the quantum ground state is described by the exponential law, poo ~ e~'"^* 
with r the decay rate per unit volume and per unit time. It consists of contribution 
from each band r„; F = X^^i r„. Precise formula from each band is given in Ref 



13 1, and goes like 



oc m|, 



For clarity we included the mass m dependence of field. This formula implies that 
the process from the n— th band is a multi-particle reaction of — * 2^9 with (y9 
particle energy of nm^/2. The most obvious contribution is from the first band, for 
which Fi = (yf^m| / (647r) = 7^ with 7^ one particle decay rate of the infiaton 
(fiip) and = ^^/2 the infiaton number density. This relation confirms the 
expected result that the leading effect in the small ^ amplitude limit is the stochastic 
2-body decay. 

Particle production accompanies back reaction against the infiaton and ultimately 
the infiaton oscillation ends by damping due to particle production. For a global 
consideration this back reaction can be treated via energy balance between the 
infiaton and the created particle energy. With cosmological expansion included, time 
evolution equation is given by 

^ + 3//p, = -|iV(p,)-7,p,, 
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= yG(p£ + p,), 

Here is a number of bose fields contributing to the parametric resonance, 8 in 
the minimal supersymmetric model in which 2 Higgs doublets couple to the infia- 
ton field. We also considered fluctuation {ip^) associated with particle production 
(fluctuation-dissipation theorem), causing the ^ mass variation. It turns out that 
the source term on the right hand side {p^p) can be handled without invoking 
integro-differential equations both in the large and the small amplitude limit. The 
remaining complexity is how to match these two regions. We used two types of 
matching, one crude and the other somewhat elaborate. Result is not very sensitive 
to these details of the matching. To understand essential features of our result, it is 
sufficient to consider the simpler matching, based on the quartic inflaton dominance 
until 6y = 2g^/m^ = ^/152 — 12 ~ 0.33 and the subsequent Yukawa dominance. 

It is now appropriate to discuss consequences that come out of numerical inte- 
gration of evolution equation. We mainly studied this under the initial ^ dominance; 
P((to) = I "^1^0' Pr(to) = 0. Evolution equation was integrated until almost com- 
plete Born decay, conveniently taken at t = = 10 7^^ . What typically occurs is 
that after O[10 — 100] x m^^ oscillations catastrophic particle production terminates 
the initial phase of dissipative ^ evolution and subsequently enters into a much slower 
phase under the small amplitude oscillation, ultimately ending with the Born decay. 
It can be checked that strong interaction leading to thermalization occurs immedi- 
ately after the catastrophic particle production. The effect of parametric resonance 
in the large amplitude regime is found most significant for a range of couplings, 
g = 10"'^ — 0.1 . In this range of coupling the initial thermal temperature denoted by 
Tj is much higher than what the Born rate would suggest; Tj = O[10^^ — 10^^] GeV. 
This high temperature phase is followed by the much slower phase of the Born decay, 
until the final temperature becomes of order 10^^ g GeV , at the time te- This agrees, 
remarkably well, with the naive estimate assuming instantaneous decay and instan- 
taneous thermalization due to the Born decay, Tb ~ 0.1 ^T^m^ ^ 10~^ (7 y^m^m^. 
But it should be noted that this simple estimate completely ignores the parametric 
resonance effect, which predicts an epoch of much higher temperatures: the thermal 
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history after inflation cannot be described by a single reheat temperature such as 
Tb. 

Existence of the high temperature phase necessitates to reconsider high energy 
processes such as baryogenesis and gravitino production. Indeed, what we found in 
is that high enough initial temperature favored by GUT baryogenesis is realized 
even with low enough gravitino abundance consistent with nucleosynthesis. The old 
formula of the gravitino abundance ~ 10~^ ^ with s the entropy density is 

not valid since the thermal history right after the catastrophic particle production 
is complicated in some region of the coupling g and Tb is not a good measure to 
characterize the thermal history. Let us then consider gravitino production after 
inflation, following evolution equation for the gravitino number density n3/2, 

where is the thermal number density of one species of created particles. The cross 
section (St>) of gravitino production ip(p ^ 93/293/2 is different 0, whether or not the 
gravitino is lightest supersymmetric particle (LSP). For the heavy gravitino (St") ~ 
250/mpi. Both possibilities of stable and unstable gravitino remain viable. Let us 
only mention a possibility of the gravitino dominated universe at the present epoch. 
With the initial temperature Tj > 2 x 10^^ GeV imposed to give a favorable situation 
for GUT baryogenesis, there is a region of parameters for the closure density of 
gravitino dominated universe if 7713/2 = 0.1 — 10 GeV. The basic reason this becomes 
possible is that the initial large gravitino yield created right after the catastrophic 
particle production is much diluted via the late phase of Born decay. This is again 
a reflection of the large disparity of the two temperatures; Tj ^ Tg . Of course, it 
remains to demonstrate a sizable baryon to photon ratio. But things are not bad: 
there is an epoch immediately before the catastrophic particle production in which 
non-equilibrium environment necessary for baryon generation exists, and moreover 
the baryon to photon ratio is of order lO"^'' allowing some amount of dilution in 
later epochs. 
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